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ABSTRACT 



Observations suggest that relativistic particles play a fundamental role in the 
dynamics of jets emerging from active galactic nuclei as well as in their interaction 
with the intracluster medium. However, no general consensus exists concerning 
the acceleration mechanism of those high energy particles. A gravitational ac- 
celeration mechanism is here proposed, in which particles leaving precise regions 
within the ergosphere of a rotating supermassive black hole produce a highly col- 
limated flow. These particles follow unbound geodesies which are asymptotically 
parallel to the spin axis of the black hole and are characterized by the energy 
E, the Carter constant Q and zero angular momentum of the component L z . If 
environmental effects are neglected, the present model predicts at distances of 
about 140 kpc from the ergosphere the presence of electrons with energies around 
9.4 GeV. The present mechanism can also accelerate protons up to the highest 
energies observed in cosmic rays by the present experiments. 

Subject headings: Kerr geodesies, astrophysical jets, high energy particles 
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1. Introduction 

Multiwavelength observations of different astrophysical objects indicate the presence 
of "jets" (Marscher 2005) driven probably by a compact object like a black hole (BH) 
or, in some cases, a highly magnetized rotating neutron star. Jets are observed at scales 
ranging from sub-parsec up to hundreds of kiloparsec. Gamma-ray bursts are an example 
of small scale jets, since they are supposed to be the consequence of shocks occurring in 
highly collimated relativistic flows (Meszaros et al. 1999; Frail et al. 2001; Rossi et al. 
2002) originated either at the death of a massive star or when a neutron star merges with 
another neutron star or with a black hole. Large scale jets are observed in association with 
radio-galaxies, quasars, blazars or AGNs (active galactic nuclei) in general and their origin 
is probably the consequence of the twisting of magnetic fields anchored in the very inner 
region of an accreting disk around a supermassive black hole (Blandford & Levinson 1995; 
Meier et al. 2000; Camenzind 2005). 

Jets associated with AGNs have a complex structure with bulk motions characterized 
by Lorentz factors typically in the range 10-50 and total power ranging from 10 44 up to 10 47 
ergs -1 . The composition of these flows is still uncertain but, in general, supposed to be 
constituted by electrons and protons or/and electron-positron pairs or even heavy nuclei. 
The composition of the jet is certainly related to the physical processes that create and 
energize the flow, representing an important key for the understanding of the launching 
mechanism. The large scale acceleration required to explain the high velocities of the bulk 
motion cannot be purely hydrodynamic and are probably a manifestation of the presence 
of extended magnetic pressure gradients (Vlahakis & Konigl 2004). Fully relativistic 
simulations of accretion disks around a rotating black hole indicate that unbound flows 
can emerge self-consistently from the accretion flow (De Villiers et al. 2005). According 
to these simulations, the flow has two main components: a hot, fast and tenuous outflow 
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along the jet axis and a cold, slow and dense flow along the funnel wall defining the jet 
geometry (see also Meliani & Keppens 2009). For slow rotating BHs, the flow energetics 
is dominated by the kinetic energy and the enthalpy of matter whereas for fast rotating 
BH the energetics is essentially given by a Poynting flux. Jets dominated by kinetic energy 
penetrate easily in the intra-cluster medium (ICM), forming low density cavities elongated 
in the radial direction (Guo & Mathews 2011). This is not the case if the energetics of the 
jet is dominated by highly relativistic particles. In this case, due to the low inertia, the 
jet decelerates rapidly in the ICM, producing large cavities due to the lateral expansion 
produced by the pressure of the relativistic particles (Guo & Mathews 2011). 

The result of these simulations emphasize the importance of the presence of relativistic 
particles in the jet, either to characterize the dynamics of the flow or the interaction 
with the ICM. However, the acceleration mechanism (or mechanisms) of these relativistic 
particles is not yet well understood although magnetohydrodynamic shocks and Fermi-like 
mechanisms are often invoked as possibilities. Here we examine a completely different 
alternative, i.e., a purely gravitational acceleration process based on the presence of an 
ergosphere around a Kerr BH. In our scenario, we assume that an accreting rotating BH is 
present in the center of an active galaxy. Matter penetrating the ergosphere can undergo 
the Penrose process (Penrose 1969) and, under certain conditions, the emerging particles 
follow geodesies asymptotically parallel to the rotation axis, acquiring very high energies. 
Although the efficiency of the Penrose process is still a matter of debate, this important 
question will be not examined in the present paper. 

The existence of unbound geodesies leaving the ergosphere along the z-axis and focusing 
at infinity was already demonstrated by Gariel et al. (2010, hereafter GMMS10). In the 
present work, it is shown that test particles, independently of their electric charge, following 
these highly collimated geodesies, model a narrow energetic beam in precise regions of the 
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ergosphere. The present investigation is addressed to the study of some particular solutions 
describing those geodesies as well as to the analysis of the initial conditions required for 
their existence. The paper is organized as follows: in Section 2, the constants of motion, 
the total particle energy E and the Carter constant Q are derived as functions of two real 
roots of the characteristic equation R 2 (r) = 0, where the function R governs the timelike 
geodesies in the Kerr's metric; Section 3 is dedicated to the study of the particular case in 
which a double root of the equation R 2 {r) = exists. It is shown that high particle energy 
values are possible only for two narrow ranges of values of the considered double root; in 
Section 4 the two remaining roots of the characteristic equation are examined as well as the 
consequences for the allowed values of E and for the asymptotes of the geodetic motion; in 
Section 5 an analysis of these various solutions is performed, permitting to restrict to one 
the different possibilities and, finally, in Section 6 the main conclusions are given. 



Firstly, it should be emphasized that our model is highly idealized since interactions 
with the ambient medium or with magnetic fields, which affect the motion and the energy 
budget of particles, are not included in the present approach and will be considered in a 
future paper. This investigation will be focused on the study of the motion of test particles 
following unbound geodesies along the z-axis, under the assumption that the central 
supermassive black hole (SMBH) rotates steadily. 

Assuming an axisymmetric geometry, the generalized cylindrical or Weyl coordinates 
(p, z, 0), related to Boyer-Lindquist generalized spherical coordinates (r, 9, 0) by 



2. 



The constants of motion 




(1) 
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where 



\MJ 



(2) 



are the most suitable for describing the system. As already mentioned, the existence of 
special unbound geodesies in this frame was recently demonstrated by GMMS10. These 
geodesies stem from the ergosphere and when z — > oo, they are asymptotically parallel to 
the z axis and positioned at distance 



that depends on p e = a/M and on two constants of motion; the Carter constant Q and the 
energy E, while the third constant of motion, the z component of the angular momentum 
L z is necessarily null. 

The function R(r) (see, for instance, Chandrasekhar 1983) introduced in the expression 
of the Kerr timelike geodesies (test particle mass \/5\ = 1) plays a fundamental role in 
the analysis of the jet collimation when the engine at the centre of the accretion disk is 
supposed to be a stationary rotating SMBH. This function is a fourth order polynomial, 
i.e., 



where the coefficients, excepting 03, depend on the constants of motion and on the BH 
parameters (Chandrasekhar 1983) as 



Without loss of generality, we have put M = 1 and L z = when considering the 
special 2D-geodesics given by Equation (jSJ). Hence, the spin a of the SMBH being fixed 
(— 1 < a < 1), we have two independent parameters left, Q and E, or, equivalently from 
Equation (j3J), the position p x of the asymptote parallel to the z-axis and the energy E. 




(3) 



R 2 (r) = a^r 4 + a^r 3 + c^r 2 + a\r + a 



(4) 



a = -a 2 Q, a x = 2{a 2 E 2 + Q) 
a 2 = a 2 (E 2 - 1) - Q, a 3 = 2, a 4 = E 2 - 1 



(5) 



-7- 



Let us consider the possible roots of the equation R 2 (r) = of the characteristics r = 
of the autonomous system of geodesies equations (Chandrasekhar 1983), i.e., 

a 4 r 4 + a 3 r 3 + a 2 r 2 + a x r + a = (6) 

The polynomial given by Equation (E]) has four roots, labeled with z = 1,2,3,4, which can 
be a priori real (positive or negative) or complex (contrary to the physical variable r, which 
is always real and defined in the interval 1 + \/A, oo ). The two equations R 2 {r{) = and 
R 2 { r 2) — are linear in Q and (E 2 — 1). Thus, the solution of this linear system permits to 
express the constants of motion as functions of the roots r 1 and r 2 , namely 

Q = ^p- K + a * h(n - 2) + r 2 (r 2 - 2)] + r\r\} (7) 



and 



{E 2 - 1) = ~ {a 4 + a\r\ + r 2 2 ) + r x r 2 [n(r 2 - 2) - 2r 2 ]} 



with the denominator D given by 
D = a 4 (2 + ri+r 2 )+a 2 [r 3 + r 2 r 2 + r x r 2 {r 2 - 4) + r 3 ] +r x r 2 [(r 2 + rir 2 )(r 2 - 2) - 2r 2 ] . (9) 

In the third possible equation, R 2 (r^) = 0, the parameters Q and E 2 — 1 can be 
replaced by Equations (CO) and (jHJ), leading to a relation between r 3 , ri and r 2 , allowing 
in principle, to determine the values of r 3 as a function of ri and r 2 only, since the spin 
parameter a is fixed. The fourth possible equation, -R 2 (r 4 ) = 0, does not represent any new 
result since the roots r 3 and r 4 are the same. 

In Equations (j7J) and (jHJ), it is worth noting the symmetric role of r% and r 2 , and that 
Q and E 2 — 1 have the same denominator D. Thus, if and only if, it cancels, we have both 
E — > oo and |Q| — > oo, whereas pi, since it depends only on their ratio (see Equation (j3J)), 
tends towards a finite value. From Equations (J3J), (JTJ) and (jSJ) we obtain for the position of 
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the asymptotes 

(P±) 2 = (a 2 + r 2 )(g 2 -r 1 r 2 )(g 2 + r|) 

\p e J ' a 2 {a 4 + a 2 (r 2 + r 2 ) + rir 2 [ri(r 2 -2)-2r 2 ]} 

It should be emphasized that the roots r\ and r 2 , satisfying the condition D = 0, for 

a given value of the spin parameter a, permit to define regions in the ergosphere from 

which unbound geodesies parallel to the z-axis emerge, along which move particles with 

"quasi" -infinite energies and that could eventually explain the very energetic particles 

observed in cosmic rays. We will return to this point later. 



3. Real roots n = r 2 

In order to simplify the mathematical analysis and without loss in the physical insight, 
we assume in this paper the particular case in which a double real root r\ = r 2 = Y exists. 
Under this assumption Equation (Tj0) can be recast as 

R 2 (r)=a 4 (r-Y) 2 (r 2 + Br + C) (11) 

and Equations (j7j) and (JSJ simplify as 

_ [a 4 + 2a 2 (y-2)y + r 4 ]y 2 

^ a 4 (l + Y) + 2a 2 (Y-l)Y 2 + (Y-3)Y 4 1 ' 

and 

F 2 1= a 4 + 2a 2 F 2 + (F — 4)F 3 

a 4 (l + F) + 2a 2 (F-l)F 2 + (F-3)y 4 1 J 

When S — > oo, of course we have also \Q\ — > oo but, as already mentioned, their ratio tends 
towards a finite value in order that the coordinate p% remains finite, i.e., 

>iV (a 2 -F 2 )(a 2 + F 2 ) 2 



p e y a 2 [a 4 + 2a 2 y 2 + ( Y - 4) Y 3 ] 

In order to perform some numerical estimates, we will assume, unless otherwise stated, 
a "moderate" rotation for the SMBH, fixing a = M/2 (a value adopted already by GMMS10 
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in their investigation). Hence, the functions E 2 — 1 = F(Y) and (pi/p e ) 2 = G(Y) can be 
plotted as shown in Figures 1 and 2. 

Since F(Y) and G(Y) must be simultaneously positive, the only possible solutions 
correspond to the two intervals 

Ye[-0.5,Y 0a ] (15) 

and 

Y e [Y ob , 3.8697] (16) 

with Y 0a and Yq& being asymptotes of F(Y) for which E — > oo. These asymptotes can be 
evaluated numerically since they correspond to the roots of the equation D = 0. For the 
assumed BH parameters, it results Y 0a ~ —0.2418 and Y ob ~ 2.8832. 

Hence, there are only two possible values of p\ for which E — > oo, corresponding to the 
two intervals defined above. For these limits, Y = Y^ a — e and Y = Y~ob + e i when e — > 0, we 
obtain respectively for the coordinate p\ 

— ~ 0.6932 and ^ ~ 10.2411. (17) 

Pe Pe 

For the upper bound of the interval (fTBl) . i.e., for Y = 3.8697 where p\ — > oo, we 
have E 2 — 1 = and, for the lower bound of the interval (1151) . i.e., for Y = —0.5 where 
E 2 — 1 = 2, we have p\ = 0. 

When n 7^ r 2 , the problem is quite difficult and a more careful study is necessary. This 
is presently underway and will be reported in a future paper. Nevertheless, the following 
aspects may be anticipated. As we have seen, in the case of a double root, geodesies 
followed by particles of sufficient high energy (E > 3), remain clustered around values of p\ 
corresponding to E — > oo. Moreover, in the general case the condition D(ri,r 2 ) = defines 
the regions from which high energy geodesies emerge. The later condition is given by a third 
order polynomial in r 2 , which has a unique real root that can be explicited as a function of 
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Fig. I. — Plot of E 2 — 1 = F(Y), where E is the energy of the test-particle as a function 
of the double root Y, computed from Equation (13) for a BH of mass M = 1 and spin 
parameter a/M = 0.5. The values of Y for which E 2 — 1 is positive, as expected for unbound 
geodesies, are clearly seen as well as the three values of Y for which E 2 — 1 tends (positively) 
to infinity. 
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Y 

Fig. 2. — Upper panel: Plot of the ratio (pi/p e ) 2 as function of the double root Y, evaluated 
from Equation (THl) . for a BH of mass M = 1 and spin parameter a/M = 0.5. The intervals 
of Y in which the ratio G(Y) is positive can be seen. Lower panel: details of the variation 
of the ratio G(Y) corresponding to the rectangle shown in the upper panel. 
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7*1. This can be inserted into Equation f lTUj) defining the position of the asymptote p\ as a 
function only of r\. That leads again to very restricted ranges of possible values of p\. 

For the sake of completeness and in order to investigate also the possible influence of 
the spin parameter in our analysis, we have also considered the case of a extreme Kerr black 
hole (a/M = 1). The results are qualitatively the same. The only physical solution leading 
to E — > oo corresponds to the double root F=-0.4142 and to an asymptote whose position 
is pi ~ 0.8284. 



4. Roots r 3 and r 4 

Identifying Equation (II ip with Equation (j3J), rewritten in terms of the parameters 
E 2 — 1 = F(Y) and (pi/p e ) 2 = G(Y), without the explicit form of these functions of Y (see 
Equations ( !T3l) and ( fl4|) ). yield the four relations, 

B-2Y = -, G+j = 2Y(BY-2C), 
1-G = 16CY 2 , 2-G = 4(C -2YB + Y 2 ), (18) 

which are linear in 1/F, G, B and C. After eliminating 1/F and G, we obtain from the 
equations above 

2(AY 2 -1)[Y(AY-1) + Y-1] 

(AY 2 - l) 2 - 16r 2 (4F-l) ' 1 ' 

r (4F 2 - l) 2 + 16Y(Y - 1) 

4[(4F 2 - l) 2 - 16F 2 (4F - 1)] ' 1 J 

Hence Equation fllip can be recast as 

R 2 = a A (r-Y) 2 (r 2 -Sr + P) (21) 
= a 4 (r — Y) 2 [r — r 3 )(r — r4), 

where r 3 and r 4 are the remaining roots, in general distinct, and we have introduced 

S = r 3 + ri = -B, P = r 3 r 4 = C, (22) 
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B+ (B 2 -AC) 1 ' 2 , (23) 
B- (B 2 -AC) 1/2 ] , (24) 



or 

1 

r 3 = — 
3 2 

1 

r4 = "2L 

where -B(y) and C(Y) are given respectively by Equations (|T9|) and (|20|) . 

The curves r 3 (F) and r 4 (y) are real only for some well defined ranges of Y. In 
particular, in the interval ( I15p defining Y, r 3 and r 4 are not real. In order to have the 
expression r 2 + Br + C in Equation All]) real, where B and C are real, r 3 and r 4 have 
to be complex conjugated, i.e., r 3 = z = B\ + %C\ and r 4 = z. Hence, the sign of the 
expression r 2 + Br + C = (r + _B 4 ) 2 + C 2 is always positive, and P = C = B 2 + Cj>0 and 
3=-B= -2B l < 0. 

In the interval (1161) . the two roots r 3 and r 4 are real, P = C is negative (which means 
two roots of opposite signs) and B = —T is positive. The most precise value we have 
numerically obtained for the left limit (where, in principle, E — > oo) of the range f lT6|) allows 
us to reach the value E ~ 1.1 x 10 32 . Then, the corresponding real values of the roots are 
respectively r 3 = —5.8988 and r 4 = 0.1324. 

Also, it is worth observing that E is steeply decreasing either for a weak variation e 
of Y from Y b (e > 0) or from y a (e < 0), while pi is weakly increasing for this same 
small interval of Y. For example, when Y goes from Yo& to 2.922, the energy E is steeply 
decreasing from 10 30 to 3, while the position of the asymptote p\/p e increases by a small 
amount from 10.24 to 10.68, which means a large concentration of the most energetic part 
(the "spine") of the beam immediately near at the right hand side of pi/p e = 10.24. At its 
left hand side, there is no beam produced. 

Likewise, for the interval fll5p . the energy E is very steeply decreasing from "infinity" 
to 6. Within our numerical precision, when Y — > Yo a , the highest value obtained for the 
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energy was E max ~ 2 x 10 30 . The corresponding asymptote pi/p e inside the ergosphere 
decreases very slightly from 0.6932 to 0.6764. Here the jet is still more concentrated just at 
the left of coordinate value 0.6932, while beyond its right side, there is no beam at all. 

As a result, the present model predicts a radial structure of the "jet", with a well 
defined profile for the energy (or velocity) distribution of the particles. 



5. Unbound geodesies for high energy particles 

Now let us consider the set of possible energetic geodesies framing a jet and their 
corresponding asymptotes, satisfying the conditions considered in the previous section. In 
this case, the choice of admissible initial conditions (IC) depends strongly on their positions 
relative to the characteristics of the system of geodesic equations (see Equations (2) and (3) 
in GMMS10). Indeed, each characteristic separates the plane into two regions and a given 
geodesic cannot cross the borders defined by those curves. In Boyer-Lindquist coordinates, 
the characteristics are defined by the equations 

f = 0, = 0, (25) 

each of which is equivalent to the equations, 

R(r) = 0, T{9) = 0, (26) 

respectively, where R{r) is given by Equation @ and T{9) by 

T(9) = -(b i p i + b 2 p 2 + b ) 1 / 2 , (27) 

with p = cos 9 and 

bo = ^~ 2 , h = a 2 , 6 4 = -(^)V-l). (28) 

The solutions of Equation (1261) . when they exist, are the roots and the roots 9i, which 
define respectively circles and straight lines from the origin. 
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In Weyl coordinates p and z, each characteristic Equation (25) is equivalent to the 

equations (see Equations (17) and (18) in GMMS10) 

R(a 2 - A)z Ra , . 

P= _, A , * = -^-. (29) 



and 



respectively, where 



apA A 
Ta 3 p . Taz 



(a 2 - A) A A 



(30) 



( M + i)V+(jL)V (3i) 



Each set of Equations ( 1291) and (1301) leads to 

cfe a 2 p 



(32) 



(33) 



dp (a 2 — A)z J 

and 

<fc (a 2 — A)z 
dp ot 2 p 

respectively, defining the two families of characteristics for the geodesies of type Equation 
( 123]) in which we are interested, namely, ellipses (corresponding to r = 0) and hyperbolae 
(corresponding to 6 = 0). Let us note that the product of the derivatives of Equations ( 132]) 
and (133]) . of these two characteristics is —1, indicating that they are orthogonal. 

In the one hand ellipses exist when there are solutions r = r» = constant of Equation 
( 132]) for any 6, with r« > 1 + v^4 or equivalently a = = constant ( because a = r — 1) 
with ctj > vA Then Equation ( 132]) can be integrated yielding 

where K\ is an integration constant. The comparison of Equation (13^]) . valid for any 9, 
with Equation (P), implies ifi = 1. 

On the other hand, hyperbolae exist when there are solutions of cos# = p, — p^ — 
constant of Equation (133]) for any r, with p 2 < 1. These are solutions of the equation T = 0, 
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when L z = 0, with 



1 - 



a 



A. 



(35) 



where we have defined A = {Q/ [a 2 (i? 2 — 1)] + (z/a) 2 }. There are two possible cases, 
namely p? = 1, then T = for any Q, or pf = -Q/[a 2 (E 2 - 1)] = 1 - (pi/p e ) 2 < 1, 
being positively defined only if Q <0, or equivalently, if pi < p e . Then, for p 2 = 1, we 
have z = r — 1 = a and p = for any r and Equation ( 13"3"|) reduces to — >• oo , and the 
characteristics being along the semi-axis z > v^4. For p 2 = — Q/[a 2 (£' 2 — 1)] we have from 
Equation (|33|) 



dp 

which can be integrated leading to 

P = K 2 



z 1 - M 

zp 



(36) 



1 1/2 



i-l -a 

p. 



(37) 



where is an integration constant. Comparing Equation (!37|) . valid for any r, or 
equivalently for any a, with Equation (CQ) yields K 2 + p 2 = 1. 



Relation (l3Tj) represents a family of hyperbolae parametrized by 



Pi 

Pe 



2\l/2 



yielding 



1 

A 



1- ^ 



n -1 



2 1 / Pi 1 2 



1 



(38) 



(39) 



If the IC of a geodesic lies inside an ellipse of the type defined by Equation (j34j) , such 
a geodesic cannot be unbounded and hence cannot go to infinity. So, the allowed IC has 
to satisfy a triple condition, i.e., i) be inside the ergosphere in order to be possibly issued 
from a Penrose process; ii) be outside the largest elliptic characteristic, corresponding 
to the larger value of the roots r« and iii) be above the higher hyperbolic characteristic, 
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which corresponds to the highest value of the roots \pi\ < 1. These conditions restrict the 
permitted domain of IC. 

An ellipse, solution of Equation when it exists (i.e. when G [1 + vA,oo[), 
can intersect the ergosphere only if its semi-minor axis 6, = (a 2 — A) 1 / 2 is smaller than 
p e = a/M, i.e. if < 2. 

As an example, let us consider again the special case of a double root Y studied 
previously. 

a) The first admissible range is Y G [Y ob , 3.869] (see Figure 2). These roots, 
belonging to the domain of physical definition r G [1 + \J~A~ = 1.86, oo[, correspond 
to the existence of elliptic characteristics. The smallest ellipse has as semi-minor axis 
bi = [(Y - l) 2 - A} 1 / 2 = [(1.88) 2 - 0.75] 1 / 2 = (2.78844) 1 / 2 along p, and as semi-major 
axis a,i = \fo?i = Y — 1 = 1.88 along z, obtained for the smallest value Yq& — 2.88, 
corresponding to p\j p e ~ 10.24. This ellipse contains the ergosphere, the limits of which 
being z mSLX = \J~A~ = 0.866 and p max = 1/2. Hence, it is always impossible to have IC 
simultaneously inside the ergosphere and outside any ellipse. Thus, in this case, it is 
impossible to have unbound geodesies starting from the ergosphere. 

b) For the second admissible range, we found that Y G [—0.5, Y 0a ] (see Figure 2). These 
roots do not belong to the domain of definition of the physical variable r, which means that 
there is no corresponding elliptic characteristic. The only remaining possible limitation 
depends on the position of the hyperbolic characteristics Equation ( |39l) . The hyperbola 
intersects the z -axis at the point with coordinates p = and zq = {A[l — (piPe) 2 }} 1 ^ 2 and 
tends asymptotically towards a straight line of equation p ~ ztan#i, with sin^ = p\j p e . 
The domain of possible IC is located between the z-axis, the ergosphere and the hyperbola. 
For example, for Y = -0.242, Pl /p e = 0.693, 6 1 = 21°, and z = 0.624(< y/A = 0.866). 
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As we have seen, geodesies characterized by E — > oo may exist. Indeed, 
trajectories defined by the differential Equation (21) in GMMS10, are well 
behaved when E — > oo (but with p 1 finite), as can be easily verified by factorizing 
E in the numerator and the denominator of the second member, leaving a finite 
quantity when E — > oo. The existence of such trajectories has been verified by 
the numerical solution of Equation (21) in GMMS10 for high energy values 
and the corresponding values of p\ as described in Section 3. Of course, the 
value of the energy will be always finite and fixed by the particular Penrose 
process which takes place at the origin of these particular geodesies inside the 
ergosphere. The main point resulting from the present investigation is that 
if very energetic particles are produced as a consequence of such a process, 
there are unbound geodesies that will be followed by those particles, which will 
permit their ejection in a collimated way along the spin axis of the black hole. 

In Figure 3 is plotted the geodesic which tends asymptotically towards the corresponding 
value p — pi — 0.347, for which the test particle has a very high (theoretically "infinite" , but 
for present calculations, we took value E = 10 6 and p\j p e = 0.693). This plot corresponds 
to the IC pi = 2.8 x 10 -6 and zi = 0.8521, which correspond to a point just inside the 
ergosphere at its top near the z-axis, i.e., near the event horizon. For the other limit 
(Y = -0.5), p 1 — 0,E— v 7 ^, #i = and z = \[A = 0.866. 

In order to illustrate the considered geometry, we have plotted in Figure 4 the relevant 
curves and surfaces like the ergosphere, the critical hyperbola, a possible unbound geodesic 
and its corresponding asymptote. Curves in black and with subscript a correspond to the 
case a/M = 0.5 considered above, while curves in red and with subscript b correspond to 
the extreme case a/M = 1. The geodesic for the case a/M = 1 was calculated for the 
following values of the constants: E = 10 6 and p 1 = 0.8284 (see end of Section 3) and for 
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nT 0,1- 
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Fig. 3. — Plot of the geodesic p(z) for a Kerr BH characterized by M = 1 and a spin 
parameter a/M = 0.5. The constants of motion and initial conditions are those given in the 
text. 
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the following IC, pi — 0.1 and z% = 0.1675. Notice that in the extreme case, the critical 
hyperbola, which delimitates with the ergosphere and the z-axis the region where unbound 
geodesies emerge, degenerates into a straight line starting from the origin and having a 
slope Az/Ap ~ 0.676. Moreover, the high energy beam of geodesies has a width, measured 
by the distance between the z-axis and the asymptote, which is about a factor of two larger 
than that obtained for the case a/M = 0.5. 

Let us return to the case a/M = 0.5. By modifying a little bit the preceding 
IC, we find other geodesies similar to the latter ones but diverging or, in other words, 
deviating slightly from its asymptote p\. For instance, taking as IC pi = 2.8 x 10~ 6 and 
Zi = 0.852086175 we obtain a geodesic that diverges almost rectilinearly from the turning 
point at p ~ 0.3466 and z ~ 50 x 10 3 attaining at z — 10 10 (or 1 Mpc) the abscissa p ~ 350 
(or 3 x 10~ 2 pc), which corresponds to a slight slope Ap/Az = 3.5 x 10~ 8 . This geodesic 
conserves its highly collimated character. If this geodesic is followed by an electron, its 
kinetic energy converges to the total energy E, attaining at z = 10 10 (or 1 Mpc from the 
BH) a Lorentz factor T = 5 x 10 4 , which corresponds to an energy 25 GeV. This same 
electron at z — 1.4 x 10 9 (or 140 kpc), where p ~ 50 (or 5 x 10~ 3 pc), attains a Lorentz 
factor r ~ 1.9 x 10 4 or an energy of the order 9.4 GeV. 

The local Lorentz factor along a given geodesic is defined by its usual expression 
T = (1 — v 2 1 C 2 ) -1 / 2 , where v 2 = v 2 + v 2 , v z = z/i and v p = p/i. The overdot means 
the derivative with respect to the proper time r. Gravitational effects are included 
self-consistently in the expressions defining the Kerr geodesies. In this way, T can be 
expressed as a function of the coordinates p and z. For the case previously considered 
(geodesic shown in Figure 3), the variation of the Lorentz factor along the geodesic can be 
expressed as a function of the coordinate z alone. The result is shown in Figure 5. Notice 
that in the beginning the energy is essentially under the form of potential energy (r ~ 1), 
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Fig. 4. — The ergosphere (labels la and lb), the critical hyperbola (labels Ha and lib), an 
example of a unbound geodesies (labels Ilia and Illb) and the respective asymptotes (labels 
IVa and IVb) are shown for black holes having a spin parameter a/M = 0.5 (black curves) 
or a/M — 1 (red curves). 
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while far from the ergosphere, the energy is essentially kinetic. 

6. Conclusions 

Taking the roots of the characteristic equation for unbound 2D-geodesics with L z = 
as free parameters, we have shown that the two remaining constants of motion, E and 
Q, of a test particle following geodesies that asymptotically tend to a parallel line to the 
z-axis, can be expressed as a function of the aforementioned parameters. In the particular 
case of a double root and assuming a spin parameter to be equal to a = M/2, restricted 
domains of asymptotes corresponding to high particle energies are found. This means that 
the Kerr metric can generate powerful collimated jets of high energy particles in some well 
defined regions inside the ergosphere. Indeed, for this special case, only two possible ranges 
of pi are possible, namely p x e [0.3382,0.3466] and p 1 G [5.12,5.34] for E G [VE,oo[ and 
E G ]oo, respectively (see Figure 2). 

It is worth mentioning that the present results are a strict consequence of the structure 
of the Kerr metric. The main approximation is the assumption of the existence of a double 
real root Y for the characteristic equation. Although we have considered in some more 
detail the case of a moderately spinning black hole (a/M = 0.5), the results are qualitatively 
the same for the case of a extreme Kerr black hole. We hope that by relaxing the double 
root hypothesis, maybe it would be possible to obtain thicker beams of energetic particles. 
An investigation is currently in progress and preliminary results are encouraging. 

Two positions of the asymptote corresponding to an "infinite" energy in our model are 
pi ~ 5M and p\ ~ 0.3466M. The latter is the only case compatible with the limitations 
imposed by the characteristics of the system of geodesic equations, according to our 
discussion in Section 5. The consequence of these mathematical constraints is that the 
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Fig. 5. — Lorentz factor as a function of the coordinate z for the geodesic shown in Figure 
3. 
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resulting jet is very narrow. However, there is some observational evidence for the presence 
of radial flows in some jets (Giroletti et al. 2004) and, as already mentioned, they may 
have a two-component flow, i.e., a relativistic powerful inner jet and a slower, less powerful 
outer flow (De Villiers et al. 2005; Xie et al. 2012). The present model could be related 
to the inner flow which carries most of the power, being constituted mainly by relativistic 
particles. 

Our results can also easily be extended to particles other than electrons, for example 
protons or heavy nuclei. This does not modify the "geometry" that we obtained, i.e., the 
positions p\ of the jets, but their energy only. For a proton {c 2 \f5[ ~1 GeV) the maximal 
energy we can here numerically calculate (which is theoretically as large as wanted) is about 
E ~ 5.6 x 10 25 eV = 5.6 x 10 7 EeV, which largely includes the highest energies of the current 
observed UHECR (Dermer et al. 2009; Hoover et al. 2010). Thus, our model could be 
relevant to explain not only a collimated flux of relativistic particles but also the production 
of these very energetic particles that could be related to UHECR. A detailed analysis of the 
Penrose process and of its efficiency is beyond the aim of this paper. Nevertheless, we would 
expect that the rate of emerging particles would be proportional to the accretion rate and 
to the ratio between the volume of the region in the ergosphere where unbound geodesies 
exist and the total volume of the ergosphere. A simple example of a Penrose event could 
be the ionization of the inner shells of an iron atom inside the ergosphere, with the nucleus 
being captured by the BH and the electron being ejected with a high energy. The accreted 
gas near the ergosphere has temperatures around 10 6 K (Montesinos & de Freitas Pacheco 
2011), high enough to ionize the K-L shells of iron, providing a theoretical basis for such a 
possibility. 

Recent results of the Pierre Auger Observatory (Roulet 2009) suggest a correlation 
between UHECR above 57 EeV and nearby (< 70 Mpc) AGNs. Since in the present 
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model the collimation of the relativistic particles occurs along the spin axis, it would be 
interesting to investigate if the associated AGNs are of type I or II, since in the context of 
the "unified model" these classes differ only by the inclination of the jet axis with respect 
to the line-of-sight. 

Let us recall briefly that since our model does not require magnetic fields, it can be 
applied to neutral particles as, for instance, neutrinos. If they have a mass of y/5[ = 0.33 
eV (Steidl 2009), according to our previous estimate, they could attain energies of the order 
of E ~ 2 x 10~ 2 EeV. 

Finally, different authors have recently discussed the possibility to produce high energy 
particles by collisions near a rotating BH (Grib & Pavlov 2011), raising the interest of 
having further investigations on the Penrose process (see also Bahados at al 2011). The 
present investigation can be seen as an additional contribution to this debate. 
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